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Abstract 

We numerically evaluate the expected performances of an optical frequency standard 
at 729 nm based on a single calcium ion. The frequency stability is studied through 
the Allan deviation and its dependence on the excitation method (single Rabi pulse 
or two Ramsey pulses schemes) and the laser linewidth are discussed. The minimum 
Allan deviation that can be expected is estimated to o"y(r) « 2.5 x /^/T with r 
the integration time. The frequency shifts induced by the environmental conditions 
are evaluated to minimize the uncertainty of the proposed standard by chosing the 
most suited environment for the ion. If using the odd isotope ^^Ca"*" and a vessel 
cooled to 77 K, the expected relative shift is —2 x 10~^^ with an uncertainty of 
±4 X 10~^^, mainly due to the quadrupole shift induced by the unknown static 
electric field gradient . 

Key words: optical frequency standard, Allan deviation, systematic eff'ects, 
trapped ion. 

PACS: 32.60.+i, 32.70.Jz, 32.80.Pj, 32.80.Qk 



1 Introduction 



Thanks to the recent progress made in atom and ion coohng and trapping, laser 
stabilization and high-resolution optical spectroscopy, narrow optical transi- 
tions are considered as a basis for frequency standards. At this time, two 
kinds of experiments are under study in various groups: one uses an ensemble 
of laser-cooled neutral atoms in a fountain, an optical lattice or a BEC, the 
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other one uses a single trapped laser-cooled ion (for a recent review see [1]). 
The present work is motivated by strong progress in storing, coohng and coher- 
ently manipulating single ions in Paul traps. Together with the ultra-precise 
optical frequency measurements achieved by frequency chains and frequency 
combs, these progress lead to the realization of single-ion frequency standards, 
as for Hg+[2], Sr''"[3,4], Yb+[4] and In"*" [5], and proposed for Ca"*" [6]. Our exper- 
imental project aims to build an optical frequency standard using the electric 
quadrupole transition 45*1/2 — ^ 3^5/2 of a single calcium ion at 729 nm. Among 
the frequency standard candidates, Ca"*" possesses the major advantage that 
the required radiations for cooling and exciting the clock transition can be 
produced directly by solid state or diode lasers. In addition, the existence of 
an isotope having semi-integer nuclear spin (J = 7/2) allows to eliminate the 
first order Zeeman shift, canceling a major source of line shift and broadening. 

The performances of a frequency standard are defined by the stabihty of its 

local oscillator (a laser in the optical case) and the precision achieved in the 
observation of an atomic transition. Frequency instability is due to deviations 
from a mean frequency throughout varying probe time intervals, while fre- 
quency uncertainty is caused by the atomic frequency fiuctuations induced by 
environmental conditions and by the experimental conditions for observation. 
The quality factor of optical atomic transitions can reach 10^^, which is 5 or- 
ders of magnitude higher than for microwave frequency standards and thus 
let hope better ultimate performances than the existing standards. 

The interrogation scheme used to probe the atomic transition influences the 
frequency stability of the proposed standard via the variation of the duration 
of the probe cycle. In the first part of this paper we discuss the choice of this 
scheme, by using numerical simulations to compare single-pulse spectroscopy 
with time-domain Ramsey interferometry. Systematic effects expected for a 
standard based on "^^Ca"*" may reduce its accuracy and precision, they are 
discussed in the second part. For this evaluation, we employ the specific pa- 
rameters of the Ca+-ion trap experiment in Marseille [6] as an example, but 
the discussion is kept as general as possible to remain applicable to other 
atomic species. 



2 Frequency stability 

Frequency stability is one of the major characteristics of a frequency standard. 
It can be quantified by the Allan deviation ay{T) measured for an average time 
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ground state |g> 

Fig. 1. Transitions involved in the cooling and probing of an ion for a frequency 
standard. In practice, the cooling scheme may involve several lasers because of 
a possible hyperfine structure and/or a possible decay towards other metastable 
states. In most cases, Tcooi is of the order of 20 MHz whereas Tj^ is of the order of 
1 Hz. 

where Q = fo/Af is the quality factor defined by the ratio of the clock 
frequency over its observed hnewidth, S/N the signal to noise ratio and Tc 
the cycle time required for the interrogation of the ion. 

The schemes to probe an optical transition of a single ion consist of a prepa- 
ration stage, an excitation of the clock transition, and a final detection stage. 
The cycle time is the sum of the corresponding time durations Tprep, Tg^c, 
and Tdet- During the preparation stage, the ion is laser cooled and optically 
pumped into the internal state chosen to be the ground state (see figure 1). 
The light scattered by the cooling transition is also used for detection of the 
ion. The probing of the clock transition by the local oscillator (a laser) is done 
by direct laser excitation. During this stage, all the preparation lasers are shut 
off. After excitation of the clock transition, the cooling lasers are switched on 
again. The quantum jump method [7] allows then to know if the atom is ex- 
cited or not: the absence of fiuorescence during the detection stage proves 
that the ion is "shelved" in the metastable upper state whereas the presence 
of fluorescence signal means that the ion is still in the cooling cycle. Repeti- 
tion of this measurement as a function of the clock laser frequency allows to 
measure the transition probability distribution. In the following, we discuss 
the laser characteristics and the maximum duration required for these various 
stages. The choice of the suited interrogation scheme is essential to minimize 
the Allan deviation. 

2.1 Preparation and detection 

We suppose the ion cooled to the Dopplcr limit in an RF-trap and located 
at the center of the trap, where the RF-trapping field is minimum and has 
little influence on the ion motion. Since the trap creates a quasi-harmonic 
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potential, the motion of the ion is a superposition of oscillations at different 
frequencies due to the spatial anisotropy of the trapping device. If we sup- 
pose, for simplicity's sake, that there is only one frequency of motion uJtr/'^T^ 
in the trap, the resulting atomic absorption spectrum is composed of a central 
frequency corresponding to the atomic transition ujq /2t: and sidebands sepa- 
rated by multiples of the motional frequency {ujq ±pa;tr)/27r, {p integer). The 
sidebands are resolved if the width of each band is smaller than their mutual 
separation (the strong confinement condition [8]). This can be achieved in 
miniature traps with high motional frequencies (~ 1 MHz) for all the narrow 
transitions considered as potential basis for frequency standards {Vjn ^ ^tri 
see figure 1). The intensity of each band p in the spectrum depends on the 
oscillation amplitude X of the ion in the trap like J'^(kX) [8] where k is the 
laser wavevector of the probed transition and Jp the Bessel function of order p. 
As the functions Jp{kX) have negligeable values when p > kX, the smaller is 
kX, the less sidebands are visible. Laser-cooling the ion reduces its oscillation 
amplitude X and thus the number of observable sidebands. A major step in 
the preparation of the ion is to access the Lamb-Dicke regime which is charac- 
terized by the reduction of the spectrum to few sidebands with a preponderant 
weight on the central frequency, this regime is reached if kX < 1. 

The motion of the ion can be described by the occupation rate of the vibra- 
tional quantum levels, characterised by the mean vibrational quantum number 
(n). This vibrational state can also be characterized from the classical point 

of view by an oscillation amplitude X — \tr\J'2,{n) + 1, where the length Xtj, = 

^J h/2mLJtr measures the size of the fundamental harmonic oscillator eigcnstate 
|n = 0). This length depends on the atomic mass by l/\/m and, as an example, 
is equal to 11 nm for a calcium ion with Utr — 2nxl MHz. The Lamb-Dicke 
condition kX < 1 can also be expressed by k\tr^2{n) -\- 1 < 1. The Lamb- 
Dicke parameter r] — kXfr quantifies the ability for a system {ion-|-trap} to 
reach the Lamb-Dicke regime for a given transition. It is of the order of 0.1 for 
optical transitions {e.g. 0.095 for calcium's clock transition in the trap taken 
as example). 

In most cases, the frequency of motion in the trap is of the order of 1 MHz, 
whereas the atomic dipole transition used for laser cooling has a width Vcooi 
close to 2tt X 20 MHz. On such broad transition (rcoo//2 ^ uJtr) the Doppler 
limit for laser cooling can be approximated by the one of a free atom [8] . This 
leads to a thermal population of the harmonic trap vibrational levels charac- 
terized by (n) ~ Vcooi/'^^tr — 10. The Lamb-Dicke condition r)^2{n) + 1 < 1 
is then fulfilled by the vibrational state reached by Doppler cooling ((n) ~ 10). 
This fulfillment sets the transition free of first order Doppler effect, while the 
second order Doppler effect is very small (see section 3.4). Furthermore the 
residual distribution of occupied vibrational levels still allows to drive coher- 
ent dynamics on the clock transition, as required for the interrogation schemes 
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and confirmed numerically in the following. Since the time needed to reach the 
Dopplcr cooling limit is of the order of milliseconds while the optical pumping 
is faster than the millisecond, we can estimate Tp^ep to 5 ms. 

The duration required for the detection stage depends on the fluorescence 
signal collected on the strong dipole transition. For such transitions with a 
width of ~ 20 MHz, one can expect at least 10'* counts per second (cps) over 
a stray light level of less than 100 cps. In these conditions, 10 ms-periods are 
sufficient to acquire enough signal to decide if the atom has been excited into 
the metastable state. As a consequence, 15 ms is a realistic estimation for 
the sum of the preparation and detection contributions to the cycle duration. 
To this minimum cycle duration must be added the excitation duration time 
Tf,xc- In the following subsection, we theoretically study the minimization of 
this probe time for different excitation schemes, assuming that the total cycle 
time Tc — T^xc + 15 ms. 

2.2 Choice of the excitation scheme 

2.2.1 Evaluation of the minimum Allan deviation 

The width A/ of the observed transition and its signal to noise ratio depends 
on the laser excitation scheme. The choice of a high-frequency clock transition 
(in the optical domain) allows to reach smaller Allan deviations than the ones 
obtained on frequency standards in the microwave domain. Until now, the 
narrowest optical transition linewidth has been observed on a Hg"*" ion [9] 
and has allowed to reach a relative frequency stabihty of 7 x 10~^^ over 1 s 
averaging. In this article, we discuss possible excitation schemes independently 
of the ion implied. We introduce a reduced Allan deviation a' = (Jy{T) x fQ\/T 
to quantify the expected frequency stability of the standard, keeping in mind 
that reduced Allan deviations between 1 and 10 have already been measured 
by several groups on lasers locked on atomic optical transitions. 

Let Pm be the probability for the ion to be in the metastable state once excited 
by the clock laser. The frequency of the transition, to which the clock laser 
will be locked, is deduced from the probability Pm measured several times on 
the low and high frequency sides of the transition. This method of frequency 
discrimination requires the excitation probability p„, to be around 0.5, where 
the slope of the probability distribution is the steepest and the frequency 
sensitivity is the highest. 

Several sources of noise can limit the signal to noise ratio. Among these is the 
quantum projection noise [10] which is dominant once the technical noise has 
been reduced. The laser excitation creates a linear superposition of the ground 
{\g)) and metastable states: \/T^-p^\g) + ^/Thn\m) . During the detection 
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stage, the atomic state is projected on one of these two atomic states. The 
variance of such a measurement is Pm(l ~ Pm) and causes a minimum noise 



yPmi^ — Pm) on the transition probabihty. This can be overcome by using 
squeezed states [11], which we do not consider here. 

The maximum signal to noise ratio which can then be observed is 



which is maximum at resonance {pm = 1), where the frequency descrimina- 
tion is unefficient. Thus, maximum frequency sensitivity and maximum signal 
to noise ratio are not compatible. Additionally, the finite upper-state life- 
time leads to spontaneous decay which, for long excitation time, can reduce 
the maximum excitation probability. As a consequence, finding a compromise 
between all these incompatible requirements deserve precise studies of the 
excitation scheme. 

Two excitation schemes have been experimentally tested by several groups: a 
single Rabi pulse or two temporally separated Ramsey pulses. The first one has 
been performed on Hg+ [2], In"*" [12], Sr"*" and Yb"'"[4], and the second one on 
Hg"*" [2] and Sr+ [3] . In the following we discuss the principal features of each 
method and then compare them. To quantify the relative stability allowed by 
the discussed methods, we evaluate the reduced Allan deviation by 



For a high frequency sensitivity, we assume that the laser probes the transition 
on each side of the line, on the two frequencies corresponding to an excitation 
probability which is half of the maximum probability measured for zero 
detuning {p„i can never exceed 1/2). The deduced A/ when these probabilities 
are equal is then the FWHM of the experimental linewidth. The evolution of 
the density matrix of the two levels \g) and \m) is computed numerically. The 
atomic system is defined by the metastable lifetime Tm fixed to 1 second (F^ = 
1/Tm), and it is driven by a Rabi pulsation Q. The motion of the ion is taken 
into account by a distribution of thermal vibrational levels, characterised by 
the mean vibration number (n) and defined by P{n) = ((n) /(I + ip))^ . The 
excitation probability Pm is then an incoherent weighted sum of the probability 
for each vibrational level to be excited in the metastable state. For a given laser 
intensity, the Rabi pulsation from vibrational level \n) |n) is proportional to 
(1 — rj^n) whereas it is proportional to rj^/n + 1 for a \n) ^ \n + 1) transition. 
We choose for r] the value of 0.095 calculated for a calcium ion in the trap 
described above. As the Doppler coohng leads to (n) ~ 10, the excitation 
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probability on the |n) — >■ \n) band is higher than the \n) |n ± 1) bands so 
these last ones were neglected. If Vt is the Rabi pulsation for the |n = 0) 
\n = 0) transition, Q x Lnijf) is the one for the \n) —>■ \n) transition, L„ 
being the Laguerre polynomial [13]. The linewidth of the laser spectrum 
(FWHM) is taken into account by adding a source of decoherence equal to 
this width in the operator controlling the evolution of the density matrix [14]. 
For a laser detuning and in the rotating wave approximation the density 
matrix evolves according to 

Pgg — Pmm — {pgm Pmg) ~l~ ^mpmm (4) 

Pgm = -i^LPgm + ^^^/2 {pgg - Pmm) " ^mPgm/'^ " ^Lpgm/'^ (5) 

2.2.2 Single pulse excitation 

To avoid power broadening, the narrow transition can be experimentally ob- 
served if the Rabi pulsation VL is smaller than the linewidth and the inter- 
rogation time longer than the lifetime. But because of the finite lifetime of 
the excited state, the maximum excitation probability is low and requires 
several seconds to be reached, reducing the relative stabihty even if the ob- 
served linewidth is close to the natural width. In an ideal context where the 
experiment is not limited by the laser stability, the calculations show that the 
smallest reduced Allan variance is reached with a single pulse which should 
last at least 1 s and drive the transition with Q of the order of 2 x F^. With 
today's laser stabihty a cycle time of a few secondes for a single measure seems 
not realistic. We rather consider excitation schemes with durations inferior to 
1 second, since a cycle has to be repeated several times before a signal can be 
built up to counteract on the frequency of the local oscillator. 

In figure 2, are plotted the excitation probability at half maximum , the 
full width at half maximum A/ and the reduced stability as defined by equa- 
tion (3), versus the cycle time Tc = Texc+i5 ms, assuming a single Rabi pulse. 
These curves reflect Rabi oscillations, which show maximum excitation prob- 
abihty for T^xc — (2g + l)7r/il {q integer). For T^^c — 2g7r/Q, the excitation 
probability is minimum on resonance and shows some maximum for other de- 
tunings. The computed FWHM has then no physical significance, which is not 
relevant here as a clock is never operated with this excitation duration. The 
interesting feature is the minimum of the reduced Allan variance observable 
for the shortest Texc- Two cases with different Rabi pulsation Q (Q = lOF^ 
and fl = 100F„) are compared in figure 2. For both cases, excitation probabil- 
ity and full width at half maximum are computed for an ion whose oscillatory 
motion corresponds to the Doppler cooling limit {(n) = 10) or to the funda- 
mental vibrational state ((n) = 0). The first result to mention is that the first 
minima of the reduced Allan variance are identical for these two vibrational 
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Fig. 2. Calculations of the reduced relative stability a (c), from the probability 
at half maximum (a) and the width A/ (Hz) (b) of the excitation probability 
profile, for an increasing excitation time T^.^^ {Tc = Tg^c + 0.015 s). This stabil- 
ity is calculated for a mctastable lifetime Tm = l/r^ = 1 s and a single Rabi 
pulse ( = lOr m or $7 — 100r„j ) for two thermal vibrational states characterised 
by {n) = and (n) = 10. Dotted line: O = lOF^ and (n) = 0, dash-dotted line: 
n = WFm and (n) = 10, broken line: Q = lOOF^ and (n) = 0, solid line: 9. = lOOF^ 
and (n) = 10. On (b) and (c), the dotted and dash-dotted line are almost super- 
imposed and only one curve is plotted on the graph. 



states, for the chosen Rabi pulsation. It confirms that Doppler coohng is suf- 
ficient for state preparation. The results shown in figure 2 suggest that the 
discrepancy between the (n) = 10 and the (n) = vibrational state increases 
with the pulse duration. We have checked that in the case of the short pulses 
we consider in the following, the results are nearly the same for these two 
vibrational distributions and thus, from now on, only the cases concerning 
(n) — are dealt with. 

To illustrate the influence of the strength of the Rabi pulsation in figure 2, 
the excitation is driven by ^2 = 10 x and Q = 100 x F^. In the first case, 
the minimum reduced Allan deviation is close to 1 but requires an excitation 
of more than 300 ms. In the second case, the first minimum is reached for a 
cycle duration of 50 ms, but this shortening of the cycle duration is paid by 
an increase of a equal to 4.5. This trend is general and a further increase of 
the Rabi pulsation leads to an increase of the minimum Allan deviation as 
well as a decrease of the required cycle time. 
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Fig. 3. (a) : Minimum reduced Allan deviation and (b) : cycle duration Tc required to 
reach this minimum versus the Rabi pulsation. Computation done for a metastable 
lifetime = 1 /^m = 1 s and a single Rabi pulse by a laser of width F l (FWHM) , 
with an ion in the vibrational ground state ((n) = 0). Broken line: Tl/2'k = 0.2 Hz 
and solid line: Tl/2-k = 20 Hz. 




We now take into account the effect of a finite laser linewidth on the min- 
imum reduced Allan deviation and compare the value computed for a laser 
as broad as the atomic transition {Tl/2ti = 0.2 Hz) to the one computed 
with a 20 Hz broad laser (FWHM). In this latter case, the performances of 
the clock are greatly reduced first by the reduction of the excitation prob- 
abihty and second by the broadening of the observed transition. The first 
drawback can be overcome by the increase of the Rabi pulsation but this is 
paid by a further increase of the transition broadening. As a consequence, for 
a given laser linewidth and a given metastable lifetime, there is an optimal 
Rabi pulsation which results in a minimum reduced Allan deviation. This is 
illustrated in figure 3 where for r^/27r = 20 Hz, a' is minimum ( a'—12A) for 
100 X i$ ^ 125 X and a cycle time of 40 ms, whereas for rj;^/27r = 0.2 
Hz, (7 is minimum for Q ~ and is then equal to 1.7 but for a cycle time of 
625 ms. In the case of the broadest laser linewidth, these results confirm the 
intuitive idea that for optimum stability, the pulse length T^^c is limited by the 
laser linewidth T^^J^l — ti" and that the Rabi pulsation is then approximately 
set by the resonant 7r-pulse condition flTg^c — tt. When the laser linewidth is 
comparable to the transition natural width (r/,/27r = 0.2 s~^), the optimum 
stability is reached for a shorter excitation time Texc^Li^m) — 0.6 and a Rabi 
pulsation under the resonant 7r-pulse condition flT^xc — 2.5, which could not 
be dedeuced from the intuitive concept. 
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2.2.3 Comparison with Ramsey interferometry 

The introduction of the separated fields method or Ramsey interferometry 
[15] was soon followed by breakthroughs in high resolution spectroscopy and 
is expected to overcome the limitations met with single pulse excitation. With 
this method the line profile is recorded after two pulses of duration T such as 
QT = n/2, separated by a free evolution time T/ree- When the laser detuning 
is scanned, the profile shows Ramsey fringes resulting from an interference 
pattern and for short enough pulse duration (or a high enough Rabi pulsation) , 
the width of the central fringe is equal to l/2T^ree and is then independent 
from the Rabi pulsation. For a chosen pulsation Q, the evolution of Pm, A/ 
and a docs not show oscillations with Tc , like for a single Rabi pulse, a takes 
very high values for short Tc and decreases toward a limit for longer T^. This 
limit depends on the choice of the pulsation fi. 

When the width of the laser is taken into account by the relaxation it causes 
on the coherence, for a given Rabi pulsation, the reduced Allan deviation 
decreases for increasing free evolution time until it reaches critical time where 
the width of the laser broadens the line. This behaviour results in a minimum of 
the reduced Allan deviation reached for this critical cycle time and depending 
on the laser linewidth and the Rabi pulsation. For increasing Rabi pulsation 
Q, this minimum Allan deviation decreases towards a limit which depends 
very few on the laser linewidth as it varies by less than a factor of 2 over 
the whole range of the considered linewidth ( 0.1 < r^/27r < 100 Hz). This 
is made possible by the short interaction time with the laser, allowed by a 
strong Rabi pulsation Q,. In the case where the experiment is not limited 
by the available laser power, Q must be chosen to reach a minimum Allan 
deviation close to the limit but also to cause a negligcable light-shift on the 
atomic levels. This light-shift is evaluated in section 3.3 and our calculations 
and numerical computations show that a Rabi pulsation Q = 1000 s~^ allows 
to reach the Allan deviation limit to better than 1% and to cause a negligeable 
light-shift. Figure 4 shows this minimum reduced Allan deviation for a laser 
linewidth from 0.1 Hz to 100 Hz (FWHM), compared with the one reached 
with a Rabi excitation scheme, where the best fl and are found to reach 
the minimum Allan deviation, like explained in figure 3. The minimum Allan 
deviation expected for a narrow laser (0.1 Hz) for a Ramsey excitation scheme 
is 1.06 and for a Rabi scheme is 1.76. These values are very close to each-other 
but their evolution with the width of the laser is very different in the two 
cases. As shown in figure 4, for a Rabi excitation scheme, the minimum Allan 
deviation increases with the laser hnewidth to reach a' = 55 for ri/27r = 100 
Hz whereas it reaches 1.48 for a Ramsey scheme with the same laser linewidth. 
This evolution illustrates the great advantage of Ramsey two separated pulses 
excitation over a Rabi single pulse. 

To give an insight of how such performances are reached, the excitation proba- 
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Fig. 4. Minimum reduced Allan deviation versus the laser linewidth Tl/2'k (Hz), 
for a mctastabic lifetime = l/r^ = 1 s. The dashed line shows the minimum 
achievable reduced deviation with a single Rabi pulse (the cycle time and Rabi 
pulsation are computed to minimize a ) . The solid line shows the minimum devi- 
ation achievable with a two Ramsey pulses excitation for a fixed Rabi pulsation of 
lOOOr^ =1000 s-i . 
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Fig. 5. Excitation profile leading to minimum Allan deviation for a laser linewidth 
= 0.1 (Hz), for a metastable lifetime = 1 s. Solid line : profile reached for 
two Ramsey pulses, with $7 = lOOOrm and Tc = 313 ms. Dashed and dot-dashed 
line : profiles reached for a single Rabi pulse, two configurations lead to the same 
a = 1.76: n = 2r^ and Tc = 1.093 s (dashed line) and n = IF^ and = 392 ms 
(dot-dashed line). 

bility profiles calculated for the conditions giving the minimum Allan deviation 
for a 0.1 Hz laser linewidth are shown in figure 5. For a Ramsey scheme and 
the chosen Rabi pulsation of 1000 s~^, the optimum cycle time is 313 ms and 
results in a 1.76 Hz wide profile (A/) with = 0.465. For a single pulse, the 
minimum Allan deviation can be reached by two different Rabi pulses. One 
with fl = 2 s^^ and lasting 1.093 s gives rise to a narrow (A/ = 0.86 Hz) but 
few excited profile {pm — 0.22). Another one with Q = 7 s~^ lasting 392 ms 
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Table 1 

Minimum Allan deviation expected for different optical frequency standard candi- 
dates, computed with a very narrow laser {Tl/Itt = 0.1 Hz) and two Ramsey pulses 
with n = 1000 s-^ 



atom 


Tin (s) 


A (um) 


^ III /III 




In+ 


0.2 


236 


1.78 


1.4-iVv^ 


Ca+ 


1.15 


729 


1.05 


2.5-iVx/^ 


Hg+ 


0.08 


282 


2.78 


2.6-iV\/^ 


Sr+ 


0.4 


674 


1.40 




Yb+ 


0.05 


436 


4.04 





leads to a broader profile (A/ = 2.28 Hz) with higher excitation {pm = 0.39). 
These three profiles illustrate the compromise required between high excita- 
tion probability and narrow linewidth to reduce the Allan deviation. They 
show that different conditions can reach this compromise. In any case, the 
Ramsey excitation scheme results in Allan deviation smaller than for a Rabi 
scheme. Furthermore, this method has the advantage of keeping low deviation 
even for laser linewidth broader than the transition. 



2.2.4 Influence of the metastable level lifetime 

We have computed the minimum reduced Allan deviation for different metastable 
level lifetimes r^, for a very narrow laser {T 1/271 = 0.1 Hz) and excitation 
by two Ramsey pulses with Q = 1000 s^^ as this value allows to reach the 
limit of the Allan deviation. The results are summarized in table 1. The com- 
puted Allan deviations are very close for all the ion optical frequency standard 
candidates (between 1.4 and 5.9~^^/v^) confirming quantitatively the predic- 
tion that optical frequency standard will overtake the performances of existing 
microwave standards in the long run. 

2.2.5 Summary 

First, our studies confirm that a Ramsey excitation scheme is more appropriate 
than a Rabi one to take full advantage of very narrow atomic transition in 
the goal of building a frequency standard. They also show that the finite laser 
linewidth implies an optimum cycle time for a given Rabi pulsation, which 
can not be deduced intuitively from this linewidth as it ranges from Tc = 120 
ms for the shortest metastable lifetime listed on table 1 (r„j = 0.05 s) to 312 
ms for the longest metastable level lifetime = 1.15 s (Ca+). The Allan 
deviation expected for a Ca"*" standard is 2.5~^^/y^ which ranks well among 
the other candidates for optical frequency standard. 
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3 Frequency standard accuracy and precision 



Besides the frequency stability, the other relevant parameters defining the 
quality of a frequency standard are its accuracy and its precision. The standard 
frequency may be shifted from the atomic resonance value by any interaction 
of the atom with external fields. If this shift is constant, it only reduces the 
standard accuracy but not its precision. If this shift varies in time or can not 
be evaluated exactly, the precision is reduced also. As these effects contribute 
to the uncertainty of the future frequency standard, all the interactions of the 
ion with its surrounding must be controlled to minimise and/or to maintain 
any shift of the clock frequency. We evaluate these shifts for a calcium ion in 
order to choose the best internal state and prepare an environment for which 
these shifts are minimum. 



The ground state l^f) of the calcium ion is 

\m) is 
natural 



5*1/2) and the metastable state 



with a measured lifetime of 1152 ± 23 ms [16] which leads to a 
width for the clock transition of 138 ± 3 mHz. We require that during 
the excitation of the transition by the clock laser, all other lasers are shut off. 
This assures that there are no hght-shifts of the levels S1/2 and 0^/2 caused 
by the cooling lasers. The other major effects that can shift the standard 
frequency are due to the local magnetic and electric fields and to the intensity 
of the clock laser itself. In the Lamb-Dicke regime, the Doppler effect shifts 
the line only by its second order contribution. We first focus on the Zeeman 
effect as it governs the choice of the isotope and atomic sublevels used for the 
standard. 



3.1 Zeeman effect 



To avoid any uncontrolled or time- varying shifts, the frequency standard must 
be made as independent as possible of environmental conditions. The first or- 
der Zeeman effect can be eliminated by the use of atomic Zeeman sublevels 
with no projection of the total moment on the magnetic field. This can be re- 
alised by the use of an isotope with a half integer nuclear spin. The most 
abundant one (0.135% in a natural sample) is ^^Ca"*" with a nuclear spin 
7/2. The hyperfine structure of this isotope can be found on figure 6. Al- 
ternatively, the first order Zeeman effect could be ehminated on the ^'^Ca+ 
transition by cancellation between |5'i/2, m = 1/2^ — > D^j2,m — l/'l^ and 

D^/2-,'m — —1/2^. The exact cancellation of the first or- 
Dctter than one hertz requires a stability of the magnetic 
field better than 10^''/iT for at least a few seconds, which seems difficult to 
realise. As a consequence the use of the odd isotope appears to be the easiest 
solution to eliminate the first order Zeeman shift from the standard frequency. 



>Si/2,m = -1/2 

der Zeeman effect to 
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4 Si,? 




Fig. 6. Hyperfine structure of the levels involved in the preparation, excitation and 
detection of '^^Ca+ (/ = 7/2) [17,18]. 
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Fig. 7. Zceman shift of the metastable hyperfine sublevels of -D5/2 versus the mag- 
netic field, for Zeeman sublevel mp = 0- 

The second order Zeeman shift depends on the choice of the hyperfine sub- 
levels. We calculate these shifts by searching the eigenvalues of the Zeeman 
Hamiltonian for the |F, mj?) states. The Zeeman shift of the fundamental 



hyperfine levels S1/2, F = 3 or F = A, mp = Oy arc at least 2 orders of magni- 
tude smaller than the shift of the metastable hyperfine levels -P" = 1, 6, 
and are not relevant for the choice of the level. Figure 7 shows the quadratic 
Zeeman shifts for the different hyperfine levels of for sublevels mp = 0. 
These curves illustrate the great variation of these shifts with the hyperfine 
sublevel. Depending on the level involved, the second order Zeeman effect can 
be as large as 98.04 Hz//iT^ for F = 1 or reduced to -9.05 Hz//xT^ for F = 6 
(see figure 7). This last level should be used as a basis to reduce the second 



order Zeeman effect of the standard. As the 



S 



1/2 



D 



5/2 



transition is 
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electric-quadrupole, the selection rules AF = 0, ±1, ±2 imply that the funda- 
mental sublevel involved in the standard should be 



To ehminate the uncertainty due to the Zeeman effect, the local magnetic field 
must be kept on the 0.1 f/T level. A controlled magnetic field is still needed 
to split all the Zccman sublevels and to be able to select the mp = 0) 
Imp — 0) transition. The two closest transitions 5*1/2,4, ±1^ — > D^/2,Q,±1^ 
are split apart by ±3.5 kHz/ fiT. So a magnetic field of 0.1 fiT allows to isolate 
the mp = ^ niF = transition and can be measured by the observation 
of these neighbouring transitions. Nevertheless, such a magnetic field may not 
be sufficient to maintain a high level of scattered light by the atomic system, 
as observed for other ions under study [4]. But it is possible to recover a high 
level signal by spinning of the laser polarisation [19,4] and we do not consider, 
at this level, this reduction of the signal as a limitation. Without any magnetic 
field applied, the local field caused by the earth and the experimental setup 
is of the order of 10~^ T and such fields can be produced by 1 A in Helmholtz 
coils. Furthermore, magnetic field fiuctuations of 0.2 ^T over one day have 
been observed in an unshielded environnement [20]. As a consequence, in a 
thermalized and shielded environnement, it is technically possible with stan- 
dard current supplies of 1 A stabilized to the mA level, to compensate for 
the already existing magnetic field and to add the desired magnetic field of 
0.1±0.05/iT. In these conditions, the frequency uncertainty due to the Zeeman 
effect is 

B = 0.1± 0.05 /xT ^ dfz = -0.09 ± 0.09 Hz (6) 



Indeed, it is possible to prepare the atomic system in the yS'1/2, F — 4, mp — Oy 
state thanks to the property of dipole transition that forbids |F, = 0) — > 
\F,mp = 0) transitions. After the cooling stage, two lasers polarised parallel to 
the magnetic field (Amp = 0) and resonant with the 6*1/2, F = A and F = 3^ - 

Pi/2, -F = 4^ transitions optically pump the system in the S1/2, F = 4,mp = 
state in few microseconds, then ready for the probe stage. Actually, the cooling 
and optical pumping stage are not so simple due to a possible decay from P1/2 
to level (see figure 6), whose lifetime is the same order of magnitude as 
-D5/2 and so requires three repumping lasers, to empty the possibly occupied 
F = 3,4,5 levels. If the repumping lasers arc a{= a+ + a^) polarised, the 
cooling and optical pumping remains efficient, as long as the three repump- 
ing lasers' detunings are different from the two cooling ones, to avoid dark 
resonances [21] and as long as their polarization is spun to prevent pumping 
into dark states [19]. With its very low natural abundance, the use af such an 
isotope is technically challenging, but it has been shown that photoionisation 
processes allow to create '^^Ca"'" ions even from a non-enriched calcium sample 
[22]. 
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3.2 Interaction with DC electric fields and their gradients 



The second order Stark effect shifts the standard frequency through the cou- 
phng of the levels S1/2 and D5/2 to all the other atomic levels by electric dipole 
interaction with any DC or slowly varying electric fields. These fields also shift 
the -D5/2 level by the coupling of its electric quadrupole moment to any electric 
field gradient. In a usual miniature spherical trap, the confining electric field 
has no static component, oscillates at a frequency of the order of 10 MHz and 
its shape in the center can be very well approximated by a quadrupole. In the 
exact center of the trap there should be no oscillating field but an oscillat- 
ing field gradient. However, in a real Paul trap, patch potentials deform the 
harmonic potential well created by the RF field. They separate the minimum 
potentiel point from the zero RF-field point and static bias voltages have to 
be applied in the three directions to make these two points meet again and 
reduce any static field to less than V/cm. This step is required to be able to 
cool an ion to the Doppler limit and to reach the Lamb-Dicke regime [23] and 
can lead to an increase of the static electric field gradient. Such gradient can 
realistically reach IV/mm on 1 mm (the typical diameter of a Paul-Straubel 
trap). The local electric field is then the sum of the quadrupole oscillating 
field that traps the ion, the bias static field lower than V/cm and the isotropic 
field radiated by the vessel considered as a blackbody. Since the frequencies of 
this radiated field are far below the optical resonance of Ca"*", the field can be 
taken into account by its mean-square value averaged over all the blackbody 
spectrum, whose value is given by [24] 

(^i.)-«31.9^Q^ (7) 



in (V/m)^ with T in Kelvin. At room temperature, this field overtakes the 
static bias field resulting from compensation of patch potentials. Nevertheless, 
it can be drastically reduced by cooling the vessel, a thermalization at 77 K 
sets this field below the level of V/cm, comparable to the bias static one. 

Thanks to a symmetry property of the second order Stark Hamiltonian (which 
behaves like a second order tensor), the Stark shift of 6*1/2 is independent of the 
hyperfine level and Zeeman sublevel. As a consequence, it is also independent 
of the polarisation of the electric field and behaves hke a scalar (this property 
is true for any level with J < 1). An electric field couples the ground state 
45*1/2 to all the raPi/2 and nP^i2 levels but in the fact, the sum of the oscillator 
strength on 4Pi/2 and 4P3/2 is already equal to 1 [25] and there is no point 
taking into account other couplings to n > 4 levels. The second order Stark 
shift on 45'i/2 is then easily evaluated to —9.5 mHz/(V/cm)^. 
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The Stark effect on the -D5/2 level can be split into a scalar term, independent 
on F and nip and a tensorial part, depending on these two quantum numbers 
and on the angle 9 between the electric field and the magnetic field defining the 
quantification axis. The sum of all the oscillator strengths of the transitions 
between 3D5/2 and nPs/2{n > 4), nF^/2{n > 4) and nF-ji2{n > 4) is only 
0.48 (according to the Harvard database [25]), suggesting that there are other 
couplings with levels belonging to the continuum. Then our evaluation can 
only be a rough estimation, but it gives a correct order of magnitude. We find 
-3.9 mHz/(V/cm)2 for the scalar part and +2.1 mHz/(V/cm)2 x (Scos^^ - 
l)/2 for the tensorial part and associate an uncertainty as high as the value 
itself to take into account that there are missing couplings. The total frequency 
shift due to DC-Stark effect is then 

5fs{Sy2^D5/2,F^6,mF^0) 

= 5.6(±4) + 2.1(±2) X I^^S^lI^^ mHz/(V/cm)2 (8) 

At room temperature, the DC Stark shift is mainly due to the isotropic radi- 
ated field and is therefore 

Sfs(Si/2 ^ A/2, F = 6, = 0) = 0.39(±0.27) Hz (9) 

If the vessel is cooled to 77 K, the contribution of the radiated field is of the 
same order as the bias static field so its direction is unknown and its amplitude 
of the order of 1 V/cm. Such a field induces an uncertainty on the frequency 
of 

6fs{77K){Sy2 ^ D^/2,F = 6,mF = 0) < 12 mHz (10) 

As for the coupling of the electric quadrupole moment of the 3-D5/2 state to any 
electric field gradient, it depends on the hyperfine level, its Zeeman sublevel 
and on the shape as well as on the symmetry axis of the electric potential [26] . 
The coupling strength is due to a non spherical repartition of the electronic 
charge density and depends on the atomic orbitals of the considered level. The 
quadrupole moment ©(3D5/2) of the fine structure state can be defined as [26]: 

e(3L'5/2) = -| (3L'5/2, mj = ^\3z^ - r^\3D^/2, mj = ^ (11) 

This is calculated by considering the electronic orbital of 3D5/2 as pure 3(i 
without any mixing with other electronic orbitals. For a single electron atom 
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[27] 



2\ + 2 



(12) 



2e / 2\ 



In our case: 

e(3D,/,) = ~- (13) 

In [26], the Cowan code is used to compute (r^)^^ for Hg+. A good enough 
and simple estimation of (r^) in alkah hke ion is provided by the quantum 
defect method [28] which gives a simple relation between the energy Eni of 
the electronic level and an effective quantum number n*{l) by 

E^i^--^^a.u.,Z^Z-N + l (14) 

3^ can then be calculated using the one-electron orbital properties, with Z 
and n* instead of Z and n. This method gives for Ca"*" 

{r%, - e.6al (15) 

where Oq is the Bohr radius. The energy shift of the hyperfine sublevel |F = 6,mF 
of 3D5/2 is 



where H is a geometrical factor equal to (3cos^/3 — 1) if the field has a 
quadrupole symmetry {V (x + — 2z^), [3 being the angle between its 
symmetry axis and the magnetic field defining the quantization axis [26] . The 
frequency shift of the standard transition under investigation is then 



The hyperfine level has little infiuence on this shift as, for exemple, for the 
level F — 2, 7/11 is replaced by 17/35. With the expected gradient of 1 
V/mm over 1 mm, the uncertainty induced by this effect reaches the hertz 
level, which is high compared to the width of the clock transition in Ca"*". 
Any modification of the patch potential, due for exemple to the ion creation 
process, alters this shift and reduces the reproductibility of the standard. Still, 
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this effect can be eliminated by averaging the transition frequency measured 
with the magnetic field along three perpendicular directions, as the geometrical 
n factor is then averaged to zero [26]. The remaining uncertainty will then 
depend on the precision of the angle setting between the three measurements. 
This precision depends a lot on the vessel design and experimental setup, and 
it seems difficult to estimate this uncertainty as long as we have not performed 
the experiment. Nevertheless, other authors [29] have projected to reduce by 
50 the uncertainty induced by this shift and we assume that a reduction by a 
factor of 10 is readily achievable, which sets the uncertainty induced by the 
quadrupole effect to ±0.1 Hz. At this point, it is important to mention that 
in spherical traps, the field gradient is inferior to the one in linear traps, due 
to the confining geometry. As a consequence, in order to minimize the shift 
induced by the gradient a spherical trap is preferred to a linear trap. 



3.3 Interaction with AC electric fields 



Si/2, 4, o; 

ht- shift caused 



D5/2, 6, Oy transition is 0.75 //W/mm^, which leads to 
^J dipole coupling equal to 0.08 niHz, which is however 



During the excitation of the clock transition, only one laser is applied. It can 
still cause an AC-Stark shift (or hght-shift) of 5*1/2 and by couphng them 
to Pi/2 and P3/2 by dipole interaction or by coupling them to other Zeeman 
sublevels of and S1/2 by quadrupole interaction (the coupling with D3/2 
is far less strong). The first two couplings produce a shift proportional to the 
laser intensity 1729 equal to 1.1 x 10"*^ x 1729 Hz. The laser intensity required 
to produce the highest Rabi pulsation of 1000 considered in part 2.2.3 on 
the 
a \ig 

negligeable compared to the natural width of the transition. 

Light-shifts of a few kHz due to quadrupole interaction with other Zeeman 
sublevels have been measured on ^°Ca+ isotope [30]. In these experiments Rabi 
pulsations of 1 MHz were used with laser detunings of the order of 1 MHz. Here 
we calculate this shift in the context of the clock transition excitation for a 
Rabi pulsation equal to 1000 s~^ and a magnetic field of 0.1 /xT. The frequency 
detuning required to probe the clock transition depends on the laser linewidth 
and on the Rabi pulsation used, and is of the order of 1 Hz. By chosing ±10 
Hz for this detuning the light-shift is then not underestimated. We find an 
effect equal to ±6 mHz decreasing to ±0.06 mHz if the magnetic field is 1 /iT. 
The sign depends on the sign of the detuning. Two reasons make this effect 
very small: the small Rabi pulsation considered for such experiments and the 
small detuning required to probe the two sides of the transition (of the order 
of a few Hz). With such small detuning, the couplings of 5'i/2,4,0^ with 

D5/2, 6, 2^ and with D5/2, 6, —2^ compensate each other (and vice-versa for 

D5/2, 6, 0^ with Si/2, 4, ±2^). Nevertheless, with the laser power and magnetic 
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field values planned for the optical clock realisation, this effect overtakes the 
ones induced by dipole couplings. 



3.4 Second order Doppler shift 



The second order Doppler effect shifts the frequency transition by 
fo 2c2 



(18) 



With an oscillating ion, cooled to the Doppler limit, the velocity of the ion 
can be written like v — Vocos{u;trt + 0) and {v'^) — Vq/2. This mean-square 
velocity is calculated by Vq — uJtrX and X — \tr^2{n) + 1 (c/2.1). With the 
values chosen in 2.1, the velocity amplitude Vq is equal to 0.32 m/s leading to 
a second order Doppler relative shift given by 

^ = -2.8 X 10-19 (19) 
/o 



In the case of the Ca"*" clock transition (/o = 4.11 x 10^^ Hz) , the absolute 
shift is 0.11 mHz, which is negligeable in the reduction of the clock precision. 
This calculation confirms that Doppler laser cooling is sufficient also to reduce 
the second order Doppler effect to negligeable values. 



3.5 Uncertainty budget 



Table 2 gives the uncertainty budget expected for an atomic clock based on 
^^Ca+. At room temperature, and with the considered magnetic field, the 
major source of frequency shift and uncertainty is the Stark effect induced by 
the radiated electromagnetic field. This effect is drastically reduced in a vessel 
cooled to 77 K and then the major source of uncertainty becomes the coupling 
with the field gradient through the quadrupole moment of ,05/2 which limits 
the ultimate precision of the clock. It can be compensated by measuring the 
frequency with three perpendicular directions of magnetic field. Nevertheless, 
the obtained precision will depend on the design of the experimental setup 
and the ability to control the directions of the laser propagation and magnetic 
field. The projections made for all these major systematic shifts show that 
an atomic frequency standard based on 81/2, 4, 0^ 6, 0^ of "^^Ca"'" 

can reach an uncertainty of 4 x 10^^^, with room for improvement by better 
compensation of the quadrupole shift and better stabilization of the magnetic 
field. 
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Table 2 



Uncertaintv budget for the freqiiencv transition of \Si/-:,, 4. O) [-C5/2' 6, O) in 
^^^Ca+ ' Z ' " 



effect 


fields / conditions 


shift (Hz)@ 300 K 


@ 77 K 


second order Zeeman effect 
Stark effect 
coupled to the field gradient 
AC Stark effect @ 729 nm 
second order Doppler effect 


0.1 fiT 

radiated and bias static field 
1 V/mm^ 
0.75 /xW/mm2, 0.1 /xT 
ion cooled to the Doppler limit 


-0.09 ± 0.09 
+0.39 ±0.27 
±0.1 
±0.006 
-1 X 10-4 


-0.09 + 0.09 
< 0.012 

±0.1 
±0.006 
-1 X 10-4 


global shift and uncertaintity 




+0.3 ±0.4 


-0.09 ± 0.19 


relative shift and uncertaintity 




+7(±9) X 10-1^ 


-2 (± 4) X 10-1^ 



4 Conclusion 



We have presented a theoretical evaluation of the ultimate performances that 
can be expected from an optical frequency standard based on an electric 
quadnipolc transition of a trapped single ^^Ca"*" ion. We studied its stabil- 
ity through its Allan deviation, assuming that the signal to noise ratio would 
be limited by the quantum projection noise. Our results show that a frequency 
instability of ^ 2.5 x 10~^^ / y/r ca,n be expected. We also show that a Ramsey 
excitation scheme allows to take advantage of a very narrow transition, even 
with a laser broader than this transition, whereas this is not possible with 
a single Rabi pulse. The minimum Allan deviation is also computed for the 
other ions which are candidates for an optical frequency standard and Calcium 
ranks well within this list. 



In a second time, all the systematic frequency shifts have been estimated 
and the environmental conditions studied in order to minimize the frequency 
uncertainty. This minimization is limited by the precision reached in the suc- 
cessive orientation of 3 mutually perpendicular magnetic fields to compensate 
the coupling of the quadrupole with a field gradient. A technical chal- 
lenge for the future optical frequency standard will be to point these three 
perpendicular magnetic fields and to reduce the field gradient. In this context, 
a miniature spherical trap is more appropriate than a linear one to a frequency 
standard. Our projections show that with a first step alignement and a cooled 
vessel, a standard based on 6*1/2, 4, 0^ -D5/2, 6, 0^ of ^^Ca^ can reach an 
uncertainty of 4 x 10^^^, an order of magnitude smaller than the most precise 
actual microwave frequency standard [31]. 
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